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SUMMARY 


The National Aeronautics and Space Administration and the Federal Avi- 
ation Administration have developed concepts which save fuel while preser- 
ving airport capacity by combining time-based metering with profile descent 
procedures. NASA conducted experiments where an aircraft descended from 
cruise altitude in a fuel-efficient manner by flying with idle thrust and a 
clean (low drag) configuration. This report describes a computer algorithm 
developed to provide the flight crew with the information needed to fly 
from an entry fix (about 100 n.mi. from the airport) to a metering fix 
(about 25 n.mi. from the airport) and arrive there at a predetermined time, 
altitude, and airspeed. Additional information is calculated for the 
flight from the metering fix to an aim point near the airport. The algo- 
rithm was developed for use in an air traffic simulation to model the 
dynamics of aircraft performing flight-idle profile descents. 

The flight path is divided into several descent and deceleration 
segments. Descents are performed at constant Mach numbers or calibrated 
airspeed, whereas decelerations occur at constant altitude. The time and 
distance associated with each segment are calculated from point-mass equa- 
tions of motion for a clean configuration with idle thrust (except for one 
constant speed segment where non-idle thrust must be used). An iterative 
process is used to determine the metering-fix altitude, within a prescribed 
altitude window, and another iterative process determines the Mach number/ 
calibrated airspeed combination for the descent which allows the aircraft 
to arrive at the metering fix at a prescribed time, altitude, and speed. 

The calculation for time and distance can be simplified by using approxi- 
mate relations for variation of the rate of descent and deceleration. 

Results for the B-737 aircraft show that wind and non-standard atmos- 
pheric properties have a large effect on the flight path and cannot be neg- 
lected. Uncertainty in the descent Mach number was found to have a large 
effect on the predicted flight time, whereas uncertainty in the weight was 
insignificant. A range of combinations of Mach number and calibrated 



airspeed is possible for the descent segments leading to the metering fix. 
However, only small changes in the fuel consumed were observed for this 
range of combinations. Therefore, a combination based on flexibility for 
scheduling seems preferable. Profile descents for heavier aircraft (in the 
range of the B-747) were found to be similar to those for the B-737. 


1.0 INTRODUCTION 

The increasing cost of fuel and air traffic control problems motivated 
the Federal Aviation Administration to develop a concept called local flow 
management/profile descent. This concept saves fuel and preserves airport 
capacity by combining time-based metering with profile descent procedures. 
Presently the flight crew is responsible for the altitude and speed objec- 
tive while ATC is responsible for time delivery at the metering fix. 
However, both the pilot and controller have little or no guidance to assist 
them. A computer algorithm is developed in this report to provide the 
information needed to fly from an entry fix through a profile descent to a 
metering fix and arrive there at a predetermined time, altitude, and air- 
speed. In addition, an algorithm is included to determine the altitude of 
the metering fix such that the aircraft will arrive at an aim point within 
the terminal area with a prescribed altitude, airspeed, and distance from 
the metering fix. The algorithm was developed for use in an air traffic 
simulation to model the dynamics of aircraft performing flight-idle profile 
descents. 

The computational algorithm divides the flight path from the entry fix 
to the metering fix into the five segments illustrated in Fig. 1-1. Four 
additional segments are used to simulate the flight path from the metering 
fix to the aim point. The segments are numbered in the order of their cal- 
culation which are reversed from the order they occur. Segment 5 starts at 
entry fix with cruise altitude (h^,) and Mach number (M^.) and decelerates at 
constant altitude with idle thrust until the prescribed descent Mach number 
(Mjj) is reached. Segment 4 is a constant speed and altitude path to the 
beginning of descent. In Segment 3 the aircraft descends with idle thrust 
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and a clean configuration (low drag) at an attitude which maintains con- 
stant Mjj. As the altitude decreases, the calibrated airspeed increases. 
When the calibrated airspeed reaches a prescribed value (CAS^j), Segment 3 
stops and Segment 2 begins. Still maintaining idle thrust and a clean con- 
figuration, the aircraft descends at constant CAS^j until the metering-fix 
altitude is reached. Segment 1 is a deceleration at constant altitude 

from CAS^ to the designated metering-fix calibrated airspeed (CAS^^). 
If an ideal flight path were flown, the aircraft would arrive at the meter- 
ing fix at the prescribed time. 

The region from the metering fix to the aim point begins with Segment 
4A where the aircraft continues to decelerate until the prescribed speed 
control airspeed (CAS ) is reached. Segment 3A is a descent at CAS._ from 
hj^f to the speed control altitude (hg^.). Segment 2A decelerates the air- 
craft from CASgj. to the calibrated airspeed designated for the aim point 
(CASgp). Then segment lA is a descent from hg^. to the aim point 
altitude (h^^p) at constant CASgp. If an ideal flight path were flown, 
the aircraft would arrive at the aim point altitude at a prescribed dis- 
tance from the metering fix. 

For all the segments described above, decelerations occur at constant 
altitude and descents have constant calibrated airspeed or Mach number. 
Clean configurations and flight-idle thrust are used throughout the flight 
path except for Segment 4 where the aircraft must use non-idle thrust to 
maintain constant Mach number (M^). Time-metering aircraft to fly pro- 
file descents reduces low altitude vectoring and fuel consumption, and 
helps alleviate airport noise since the aircraft fly longer at higher alti- 
tudes near the airport. 

The Boeing Commercial Airplane Company developed an elaborate computer 
program (Ref. 1) to calculate profile descents. However, it requires a 
sophisticated computer to perform the calculations. Much simpler techni- 
ques were developed by Knox and Cannon (Ref. 2) and Knox (Ref. 3) which can 
be used on smalT programmable calculators. On the other hand, some of the 
approximations, particularly the rate of descent approximations, used in 
their analyses can lead to significant errors in some cases. 
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This report describes a computer program to calculate the profile des- 
cents described above. It is simpler than that of Ref. 1, but more accu- 
rate than the techniques used in Refs. 2 and 3. Results are presented for 
the B-737 and typical heavy-class aircraft; and the effects of winds, non- 
standard atmospheres, and optimal paths for minimum fuel consumption are 
discussed. 


5 
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2.0 SYMBOLS 


a 

^0 

^ 0 ’ ^ 1 ’ * • • > ^4 
CAS 

Cl 

0 

F 

g 

h 

L 

M 

P 

PO 

R 

S 

t 

T 

V 


Wf 

W 

X, Y 

o o 

X, Y 


a 

Y 

Y 

as 


local speed of sound, ft/sec 

standard sea-level speed of sound, kts 

coefficients for drag coefficient, Eq. 4 

calibrated airspeed, kts 

drag coefficient 

lift coefficient 

drag force, lb and distance, n.mi. 
engine thrust, lb 

gravitational acceleration, 32.174 ft/sec^ 
altitude, ft 
lift force, lb 
Mach number 

atmospheric pressure, Ib/ft^ 
non-standard sea-level pressure, Ib/ft^ 
gas constant for air, 1716 ft^/sec^/ ®R 
wing area, ft^ 
time, sec 

atmospheric temperature, ®R 
true airspeed, ft/sec 
ground speed, ft/ sec 
wind speed, ft/ sec 
fuel weight, lb 
aircraft weight, lb 

Cartesian coordinates on ground with X pointing 
North and Y pointing East, ft 

components of ground speed in X and Y directions, 
respectively, ft/sec 

temperature lapse rate, °R/ft 

inclination of flight path relative to the local 
horizontal, degrees 

ratio of specific heats (1.4 for air) 

distance interval, n.mi. 
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At 

P 

0 


time interval , sec 

3 

atmospheric density, slug/ft 

track heading measured clockwise from X-axis, degrees 
aircraft heading measured clockwise from X-axis, degrees 
wind angle measured clockwise from X-axis, degrees 


Subscri pts 

ap 

avg 

c 

d 

f 

g 

i 

max 

mf 

min 

req 

sc 

t 

TOT 

M 

X 

0 

1 


aim point 
average value 
cruise conditions 
descent conditions 
final value 
along the ground 
subscript for segment 
maximum value 

conditions at metering fix 

minimum value 

val ue requi red 

speed control segment 

value at stagnation point 

total value 

wind value 

conditions at the cross-over altitude 
sea-level or initial value 
conditions at 36,089 ft 
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3.0 ANALYSIS 


The methods used here to calculate the deceleration and descent seg- 
ments are quite different from those used by Knox and Cannon (Ref. 2). For 
these segments, the equations of motion along and normal to the flight path 
are needed. Consider the sketch below. 



Figure 3-1. Forces on an aircraft. 


Nev/ton's second_law along the flight path gives (Ref. 4): 


W ^ 
g dt 


F - D 


W sin Y 


( 1 ) 
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For level flight, y = 0 and Eq. 1 gives: 


or 


g dt 


dV 

dt 


(F 



(2) 


This equation gives the acceleration (negative of the deceleration) at con- 
stant altitude. The idle thrust (F) and drag (D) are needed to calculate 
the deceleration and they are obtained by the methods used in Ref. 1. The 
thrust is calculated from a two-dimensional interpolation of tabulated data 
of idle thrust as a function of altitude and Mach number. The drag is 
given by 

D = Cp ^ pV^ S (3) 


where p is the air density and S is the wing area. The drag coefficient, 
Cq, is obtained from fourth-degree polynomial curve fits to tabulated 
data of Cq as a function of Mach number and lift coefficient C|_. 

The coefficients Aq, A]^, A£, A 3, and A^ are determined by interpolating 
between values tabulated as a function of Mach numbers. 

For descent segments, the rate of climb is 

^ = V sin Y . (5) 

Also, the acceleration can be written as 

^ ^ dh ... 

dt " dh dt * 


Substitute for sin from Eq. 5 and ^ from Eq. 6 into Eq. 1 to 
obtain 

W^dh_p „ 
g dh dt ■ ^ " V dt 
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which can be rearranged to give 


dh _ (F - D) 
dt “ u/i . 1\ * 


( 7 ) 


This equation is used to determine the rate of climb for descent segments. 

Expressions for ^ for constant Mach and calibrated airspeed segments are 
given in Appendix C. The equation of motion normal to the flight path is 
(Ref. 4): 


W V _ 
g dt “ 


L - W COS Y . 


( 8 ) 


For small values of |y| and 


dY 

dt 


this equation can be approximated by 


L = W. 


(9) 


Since lift (L) is related to the lift coefficient (Cl) by 


L = Cl Y s, 

Eq. 9 gives the lift coefficient as 

r _ 2W 

u. - 2 • 

'■ prs 


( 10 ) 


( 11 ) 


This equation is used to calculate Cl which, in turn, is used in the cal- 
culation of the drag coefficient. 

The atmospheric density p can be obtained from the pressure and 
temperature by the ideal -gas equation of state 


P 


= _e_ 


RT 


( 12 ) 


where R is the gas constant for air and p and T are functions of altitude 
(see Appendix Aj. 
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3.1 Methods Used to Calculate Time and Distance 


Consider the segments of the profile descent without wind effects for 
now. Wind effects will be considered later. Here the total ground dis- 
tance traveled is and the track heading is e. With no wind the 

airplane heading <l> is the same as the track heading. 



Figure 3-2. Ground path for no wind. 


5 5 

Now ASjq-|- = aS^. and Atygy = ^E^ At^ for the total distance (aSjq-j-) 

and total time (At-j-Qj) required to fly from the entry fix to the metering 
fix. Four additional segments must be added to fly from the metering fix 
to the aim point. 

Consider aS^- and At^ for each segment for a prescribed value of 
^^TOT* 

3.1.1 Deceleration Segments at Constant Altitude . - The aircraft dece- 
lerates with idle thrust in a clean configuration from a prescribed CAS or 
M to a smaller value at the end of the segment. 
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From Eq. B-7 in Appendix B, the corresponding true airspeeds can be deter- 
mined from CAS or M. Then the time required for this segment is calculated 
by numerically integrating Eq. 2, i.e. 

''f 

At = /*— L = f -W_dV . ( 13 ) 

i J (dV/dt) 7 (F - D)g ^ ^ 

and the distance traveled is calculated from 


AS 




V dV 
(dV/dt) 



U V dV 
(F - D)g 


(14) 


where Vq is true airspeed at the beginning and is the true airspeed 
at the end of the segment. The integrals above are evaluated numerically 
using Simpson's rule. 


3.1.2 Descent at Constant Calibrated Airspeed . - The rate of climb is 
given by Eq. 7 as 


dh 

dt 


/I ^ 
U dh 



The right side of this equation depends on CAS and h only since is 
given by Eq. C-1 in Appendix C. Therefore, the time At^ is calculated 
numerically by using Eq. 7 in the following integral: 




dh 

(dh/dt) 





q dh 


(F - D) 



(15) 


where hQ is the altitude at the beginning and h^ is the altitude at the 
end of the segment. 
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The distance aS^- is calculated numerically also by using Eq. 7 in 
the following integral: 


AS 


i 




V dh 
(dh/dt) 



KLILiil'""’ 

(F - D) 


(16) 


3.1.3 Descent at Constant Mach Number . - The technique here is the 

same as that above except that ^ is dependent on and h, and is given by 
Eq. C-2 in Appendix C. Therefore, 


At 


3 


f — ^ = f 

J (dh/dt) J 



( 17 ) 


where h^. is the cruise altitude at the entry fix and h^ is the cross-over 
altitude which depends on and CAS^j only (see Appendix B). For this seg- 
ment V is dependent on and h only. The distance here is similarly 


AS 

3 



V dh 
(dh/dt) 



KLtli)'''"’ 

(F - D) 


(18) 


The four integrals above for the descent segments are calculated numerical- 
ly using the trapezoidal rule with increments in altitude as the indepen- 
dent variable. 


3.1.4 Constant Speed and Altitude . - This segment is actually calcu- 
lated after segment 5 because AS 5 is needed here. Since aS^q-j. is 
prescribed, then 


aS^ = AS-pQ-p - (ASp + AS2 + AS3 + ASg) 


(19) 
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and 


At 


4 



( 20 ) 


where is the velocity corresponding to M^j at entry fix altitude (h^). 

The total time for the profile descent from entry fix to metering fix 
is AtjQ-|- = At^. . For prescribed values of h^, CAS^, h^^, and 

CAS|^^, the relations above can be used to calculate ^tjQj. However, the 
proper combination of and CASd must be determined which will make 

AtjoT = At^gq. 

The specific aircraft parameters which must be known to point-mass 
calculate the aircraft's flight-idle profile-descent path are: 

1. Aircraft weight, W 

2. Wing-reference area, S 

3. Maximum and minimum Mach number and calibrated airspeed for air- 
craft in descent 

4. Idle thrust as a function of altitude and Mach number 

5. Coefficients for drag coefficient, Eq. 4, as a function of Mach 
number 

6. Fuel flow rate as a. function of altitude, Mach number, and thrust 
is needed if fuel consumption is desired. Note that Segment 4 
required non-idle thrust whereas the other segments use idle 
thrust. 


3.2 Method Used to Calculated M^/CASd Combinations for Descent 

An iterative method is used to determine the combination of Mj and 
CASq for Segments 2 and 3 which will make Atjgj = for the region 

between the entry fix and the metering fix. For a given aircraft, minimum 
and maximum values of M and CAS are prescribed. The technique used is: 


15 



(1) starting with use to calculate and 

to calculate At^^^ for this value of M^. If At^^^ < 

^^min ‘^^req ^ ^^max’ ^ P<^ofile descent cannot be performed at 
this M(j and the computer program jumps to step (4) below. If 

^^req ^ ^^max’ other means of delay would be required. 

(2) For the M^i in step (1), calculate the corresponding value of 
CAS(j which makes j At-|-g-|- - At^^^ | < 3 sec. This is done itera- 
tively using the modified regula falsi method (Ref. 5). Conver- 
gence is usually obtained in 3 iterations or less. It was found 
that this method always converged, whereas Newton's method 
diverged for some cases. The convergence criterion of 3 sec was 
chosen as a compromise betv/een the number of iterations required 
and the accuracy of the calculated results. 


At - At . 

(3) If >0.3, the current values of M. and CAS. are 

^^max ^^min ^ ° 


the combination selected. The lower limit of 0.3 means the dif- 
ference between the desired metering-fix arrival time and the 
minimum arrival time is 30% of the span of time control for the 
current value of Mj. The 30% proportionality constant was 
chosen based on the assumption that delay, rather than advance- 
ment, is prevalent. This constant could be changed to any other 
value which may be determined to be more appropriate for local 
conditions. 


(4) If >. Mj^ax’ Otherwise, increase by 0.01 and return to 

step (1). The criterion used in step (3) was chosen so that the 
profile descent would be performed at the smallest value of M(j 
which gave At^^^ to lie more than 0.3 between the maximum and 
minimum values for this M^. This allows for more delay than a 
speed up. The minimum value of generally produces minimum 
fuel consumption because the thrust required, and hence fuel con- 
sumption, in Segment 4 is minimized. 
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3.3 Effects of Wind 


It is assumed that both wind speed and direction vary linearly with 
altitude and two gradient regions may be specified. The wind affects the 
ground speed and aircraft heading as well as the time required to arrive at 
the metering fix. For the computations here, the track heading (6) is 
assumed constant throughout the profile descent. Therefore, the aircraft 
heading (tj;) will change during the profile descent. Consider the velocity 
vector diagram below 


X 



Figure 3-3. Velocity diagram on ground. 


where 

V = true airspeed 

V COSY = component of true airspeed parallel to the ground 
tj; = aircraft heading 
Vg = ground speed 
0 = track heading (constant) 

= wind speed 

= wind angle = wind direction + 180 °. 
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From the diagram, Vg is the sum of the components of V cos y and 
along the direction of Vg, 

Vg = V cos Y cos {'i^-e) + Vyj cos (0-'P^) (21) 

where (<l<-0) is the "crab" angle of the aircraft. It is assumed to be 
sufficiently small that 


cos i^-Q) 1 • 


In addition, y is sufficiently small that cos y« 1. 


Then 


V + V^^ cos (0-<l'w) 




and the ground distance traveled is 


AS = y*Vgdt = y*V dt + y*Vyj cos (0-%) dt 
= (AS)^O Ay, cos (0-'f'y,) dt. 

UTNH J 


( 22 ) 


(23) 


WIND 

This relation can be used to correct the (aS)j,q for the wind effects. 

WIND 

(Note that Vy, cos (0 - is the tail wind in the diagram.) The time 
At^' is evaluated in the same fashion as no wind except for Segment 4. 

Segment 1 

Since the altitude is constant, Vy, and 'I'y, are constant and 
evaluated at h^^^. Thus 


AS, 


WIND 


Vy, COS (0-1'yj) Atj, 


(24) 
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Segment 2 


Assume the wind speed and direction to be the same as that at the 


average altitude for this segment, i.e., h 


^mf + \ 


avg 


Then 


(25) 


WIND 


Segment 3 

^ ^ 

Here also evaluate and'f'^j at h^^g = 2 > 

AS 3 = (aS 3 )|^q + Vy cos ( 0 -^) At 3 

WIND 


(26) 


Segment 4 (Constant Altitude, h^,) 


aS^ = aSjq-j. - (aSj^ + AS 2 + AS 3 + ASg); but also 

AS = V At + V cos (0-'l» ) At . 

4 4 4 W W 4 


Therefore 



AS 

4 


V + V cos (0-<P ) 
4 W W 


(27) 


Segment 5 (Constant Altitude, h^; this segment is calculated before 
Segment 4) 

^^5 ” ^ ^^5 • (28) 

WIND 


Segments lA, 2A, 3A, and 4A are similar to Segments 1 or 2. 
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3.4 Method to Calculate Metering -Fix Altitude and Segments from 

Metering Fix to Aim Point 


The parameters input to the computer program relative to these 
calculations are: 

(1) Dgp, distance from metering fix to the aim point 

(2) Dgj., distance from metering fix to the beginning of a speed 

control segment, Segment 2A. + aS^^ 

(3) CAS at metering fix, aim point, and for Segment 3A 

(4) Altitude window at metering fix 

(5) hgp, altitude at aim point. 

The four segments from the metering fix to the aim point are shown in 
Fig. 1-1. Segment 4A is a deceleration from QAS|^f at the metering fix 
altitude to CAS^^., then the aircraft descends until the distance from the 
metering fix reaches (Segment 3A). At this point a constant altitude 
deceleration is performed until CAS^p is reached (Segment 2A), and then a 
descent to h^p. However, the metering-fix altitude (h^f) required to reach 
the aim point at the prescribed distance D^p is unknown. 

Initially, maximum and minimum values of D^p are calculated using the 
prescribed altitude window at the metering fix (maximum and minimum values 
of hjj^f). If the prescribed value of D^p lies between these bounds, the 
modified regula falsi method (Ref. 5) is used to calculate h^-j^. If the 
input value of Dgp is outside these bounds, D^p is changed to the nearest 
extremum value so that the profile descent computations for the segments 
between the entrj' fix and metering fix can be performed. 
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3.5 Influence of Parameters on AtjQj 

The separate effects of the following parameters on AtjQ^ are calcula- 
ted by changing each parameter by a small percentage and calculating the 
change in AtjQ^: 

''raf “Sd- "c *'c> ''SyOT- “• *«■ 


When linear these effects can be considered as partial derivatives of 
AtjoT respect to each of the parameters above while holding the 
remaining parameters constant. A description of the computer program is 
given in Appendix D. 


3.6 Tracking Profile Descents 

After the profile descent has been determined, and the pertinent para- 
meters printed, a tracking subprogram is called which calculates the air- 
craft's position, heading, airspeed, ground speed, rate of climb, and 
flight path angle from the horizontal for prescribed integration time in- 
tervals. The wind effects are included in this subprogram also. The 
tracking starts at the entry fix and tracks the aircraft until it arrives 
at the aim point. The accuracy of the tracking computations is dependent 
on the prescribed integration step size in time, At. 

The methods used to calculate the tracking parameters in each segment 
are the same as those described earlier to determine the profile descent. 
Additional considerations are needed, however, to calculate the aircraft 
heading with wind effects in the tracking subprogram. Since the wind 
changes the aircraft heading from the track heading, a new aircraft heading 
is calculated in each integration step. The change in heading is calcula- 
ted by the method described in Ref. 6 for a coordinated turn. This method 
involves banking'the aircraft to an appropriate angle for the turn within 
the constraints of the bank angle and time rate of change of bank angle for 

that aircraft. This part of the tracking subprogram was developed in 
Ref. 6. 
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To understand the algorithm used for the computations, the equations 


used to 
aircraft 

calculate the ground speed and track heading from the airspeed and 
heading are given below. (Refer to Fig. 3-3.) 

0 

X 

= V COSY cos 'J' + V|^ cos'l^^ 

(29) 

X 

-flit 

(30) 

o 

Y 

= V COSY sin’i' + V|^ sin'I'yj 

(31) 

Y 

= /vdt 

(32) 

e 

= tan"^ (y/x) 

(33) 

''g 

= [J2 , ?2jl/2 

(34) 

h 

■ /S « 

. (35) 

V 

■/s« 

(36) 

Eqs. 2 and 7 are used to calculate ^ and respectively, for the 

1 ast 

tv/o integrals. The integrals above are evaluated numerically using the 
trapezoidal rule. Interpolation is used at the end of a segment since time 


is the variable of integration for the tracking part. 
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At each time interval the following information is printed: 


T time from entry fix (sec) 

H altitude (ft) 

HDL track heading (degrees) 

HDG aircraft heading (degrees) 

X X-position of aircraft (n.mi.) 

Y Y-position of aircraft (n.mi.) 

VGKT ground speed (kts) 

VCALKT calibrated airspeed (kts) 

AM Mach number 

PHI bank angle (degrees) 

PHIR bank angle rate (degrees/sec) 

SS heading correction factor for wind effects (degrees) 

ROC rate of climb (ft/min) 

GAM flight path inclination angle relative to the local 

horizontal (degrees) 

3.7 Approximate Relations for Profile Descents 

The equations developed above for calculating and tracking profile 
descents require numerical integrations for the time and distance in each 
segment. For some applications these calculations may require more compu- 
tational effort than is desirable. Refs. 2 and 3 developed approximate 
relations to perform the calculations. However, their approximate expres- 
sions were found to produce significant errors in some cases. Therefore, 
different approximate relations are developed here. 
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3.7.1 Deceleration Segments . - Vq and Vf are initial and final 
true airspeeds. Assume an average acceleration, then 


Ati 



( 37 ) 


AV = Vf - Vq 


where 



(38) 


(39) 

(40) 

(41) 


The accelerations 
craft data for F and Cq. 


( 4 (4 


are calcualted from Eq. 2 using the air- 


3.7.2 Descent with Constant CAS . - Since the rate of descent changes 
slowly and monotonically with altitude, it is sufficiently accurate to use 
an average value, i.e., 



where 




are calculated from Eq. 7. 


Then 



(43) 


24 



and 


(44) 


AS. = [V 


avg 


cos (e-'^w)] At^- 


3.7.3 Descent with Constant M . - Here the rate of descent does not, 
in general, vary monotonical ly. Therefore, the average value used above is 
not accurate enough. It was found that using ^ at 3 points and performing 
the integration with Simpson's rule was much more accurate. Thus 


At. 


/ 


dh 


dh/dt 


Ah 

~ T 


/mid 


(45) 


avg 


L' ' ' ' ' J 

Properties used to calculate(^) are calculated from V 
Then ^ 


avg 


and 


ASi = [^''avg cos (e-<^w)] 4ti . 


(46) 


3.7.4 Approximate Airspeed-Mach Number Relations . - Ref. 2 gives 
approximate, yet simple, equations to relate true airspeed to Mach number 
and calibrated airspeed. They are given below. 

V(kts) = M (661 - 2.43 x 10'^ h) (47) 


V(kts) 


CAS 

1 - 0.12 X 10"^ h 


(48) 


The cross-over altitude where and CAS^j are prescribed can be de- 
termined by equating the right sides of the two equations above and solving 
for the altitude. These approximate relations are reasonably accurate for 
altitudes below 36,089 ft, and for calibrated airspeed between 210 and 
350 kts. 

The equation^ above can be differentiated to obtain dV/dh for the rate 
of climb at constant Mach number and constant calibrated airspeed. These 
relations are much simpler than Eqs. C-1 through C-4 for dV/dh given in 
Appendix C. 
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3.7.5 Approximation to Effect of Weight . - The effect of weight on 
the time can be determined from simple relations if the weight effect on 
the drag coefficient is neglected. For all segments, except Segment 4, it 
follows from Eqs. 13, 15, and 17 that 


At. 


new 


At. 


new 
old ''^old 


(i = 1, 2, 3, and 5) 


(49) 


whereas for Segment 4 


At, 


= At 


new 


AS 


TOT 


new 


4 W 
old "old 


W 


new 

''old 


(50) 


Thus the total time is approximately 


W_ 


At 


TOT 


= At 


new 


AS 


TOT 


new 




1 - 


new 

''old 


(51) 


AS 


If — « ^Hot » weight changes have only a small effect on 

4 old 

the time required to perform the profile descent. 
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4.0 RESULTS AND DISCUSSION 


4.1 Example Profile Descent 


To illustrate the application of the computer algorithm, a profile 
descent was calculated for the NASA B-737 aircraft using the following 
input data: 


Table 4-1. Input data for B-737 profile descent. 


w 

= 

87,500Jb 

VCMF 

= 

250 kts 

s 

= 

940 ft^ 

AMC 

= 

0.78 

TO 

= 

519® R . 

AMMAX 

= 

0.78 

POl 

= 

2116.2 Ib/ft^ 

AMMIN 

= 

0.60 

HC 

= 

35,000 ft 

VCMX 

= 

340 kts 

HFMAX 

= 

23,000 ft 

VCMN 

= 

250 kts 

HFMIN 

= 

19,500 ft 

DSTOT 

= 

75 n.mi . 

HAP 

= 

9,000 ft 

HDGTRK 

= 

238° 

DAP 

= 

10 n.mi. 

TREQ 

= 

700 sec 

DSC 

= 

5 n.mi. 

KSEQ 

= 

0 

VCAP 

= 

170 kts 

NO WIND 



VCSC 

= 

210 kts 





The B-737 idle thrust and drag data were taken from Ref. 1 and are shown on 
Figs. 4-1 and 4-2. Fuel flow-rate data was also taken from Ref. 1. 

A profile descent is required to take the aircraft from the entry fix 
at cruise conditions h^ = 35,000 ft and = 0.78 down to a metering fix 
which has an altitude window h„. < h . < h - and has a required 

speed of = 250 kts. The distance between the entry fix and the 

metering fix is DSTOT = 75 n.mi. and the time required to perform this por- 
tion of the profile descent is TREQ = 700 sec. 

The remaining portion of the profile descent is between the metering 
fix and the aim point. The aim point is specified to be 10 n.mi. from the 
metering fix and the speed and altitude there are CAS^p = 170 kts and 
hap = 9,000 ft. A speed control segment is specified to begin at 5 n.mi. 
from the metering fix with CAS = 210 kts. The various segments are shown 
on Figure 1-1. 

The computer algorithm calculates the segments between the metering 
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fix and the aim point first so that the altitude at the metering fix can be 
determined. From the conditions specified at the metering fix and aim 
point, it was found that the ground distance between these two positions 
must lie within the range 

43.88 n.mi. <_ DAP ^ 57.87 n.mi. 

However, the input value of DAP was 10 n.mi. Since this value is less than 
the minimum range, the computer program changes the input value of DAP to 
the minimum value of 43.88 n.mi. which gives the altitude at the metering 
fix to be the minimum value of 19,500 ft. The distance and time for each 
segment of this region are given below. 

Table 4-2. Distance and time for segments between metering fix and aim 
point. 


Segment 



lA 

2A 

3A 

4A 

Total 

(n.mi . ) 

35.54 

3.34 

1.12 

3.88 

43.88 

(sec) 

606.8 

50.4 

14.3 

45.5 

717.0 

altitude 

for segment 2A (speed 

control 

segment) was 

calcul ated 

to be 


19,150 ft. 

With the altitude at the metering fix set at 19,500 ft, the computer 
algorithm determined iteratively the combination 

M^/CAS(j = 0.67/263.7 kts 

for the descent segments between the entry fix and the metering fix. The 
cross-over altitude corresponding to these values was 

= 27,702 ft. 

The distance, time, and fuel consumed for each segment are given in 
Table 4-3. 
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Table 4-3. Distance, time and fuel consumed for segments between 
entry fix and metering fix. 





SEGMENT 



TOTAL 


1 

2 

3 

4 

5 


aS^. (n.mi.) 
At^. (sec) 

AWf . (lb) 

1.35 

25.00 

19.24 

23.13 

6.27 

74.99 

14.3 

240.6 

176.5 

215.5 

54.2 

701.1 

4.4 

72.3 

53.0 

247.1 

16.3 

393.1 


Note that the total time required is 701.1 sec rather than the input value 
of 700 sec because the convergence criterion for the iterative procedure is 
+ 3 sec. Of the 393 lb of fuel consumed, the largest amount occurred in 
Segment 4. Since this is a constant speed segment, engine thrust must be 
greater than idle thrust and, hence, the rate of fuel consumption is higher 
than for idle thrust. 


4.1.1 Influence Parameters . - For the profile descent calculated 
above between the entry fix and the metering fix, the following influence 
parameters were calculated: 

Table 4-4. Change in parameters to increase time one second. 


PARAMETER 



CASd 


^^TOT 

U 


CHANGE REQUIRED TO INCREASE 


-322 ft 
-0.0016 
-2.3 kts 
817 ft 
-0.018 
0.11 n.mi. 
25,636 lb 
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These influence parameters were used to estimate the uncertainty in At-j-g^, 
and the results are given below. 


Table 4-5. Effect of uncertainty in parameters on time. 


PARAMETERS 

%if 

f^d 

CASd 

*^c 

Me 

^^TOT 

W 


ESTIMATED UNCERTAINTY EFFECT ON At^g-p (sec) 


+200 ft 

+D.62 

+ 0.01 

+6.18 

+3 kts 

+1.30 

+500 ft 

+0.61 

+ 0.01 

+0.54 

Jp .2 n.mi. 

+ 1.86 

+5,000 lb 

+ 0.20 

WORST CASE TOTAL 

+11.31 


This table shows that the uncertainty in Mj has the largest effect (+6.18 
sec) on the total time (701 sec), whereas the weight has the least effect 
(+0.20 sec). The approximate equations given by Eqs. 49 and 50 showed that 
At, increased with weight for all segments except Segment 4 which de- 
creased with an increase in v/eight. Thus compensating effects cause the 
overall effect of weight to be small for this example. 

An investigation was also made to determine the range of linearity for 
the effect of influence parameters on the total time required for the pro- 
file descent described above. As long as AS 4 > 0, changes in AS^g^ 
affect Segment 4 only; and since the speed is constant in that segment, 
changes in AS^g^ influence At-j-g-j- linearly. Fig. 4-3 shows the influence of 
weight and cruise altitude, separately, on the total time. Changes in 
cruise altitude of + 1,000 ft produced changes in total time of less than 
1.2 sec, and the variation is nearly linear over this range. The influence 
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of weight was nearly linear for weight changes up to +10,000 lb. However, 
weight changes of j*40,000 lb (out of 87,500 lb) produced changes in total 
time of less than 1.6 sec. Fig. 4-4 illustrates the influence of descent 
calibrated airspeed and metering-fix altitude on total time. The curves 
are nearly linear for ACAS^j within ^2.5 kts and Ah^^^^ within ^500 ft. 

Changes in total time for these two parameters are greater than those in 
Fig. 4-3. Fig. 4-5 shows the influence of cruise and descent Mach numbers 
on the total time. Changes in cruise Mach number within ^0.02 produced 
changes in total time which are nearly linear. However, changes of jj3.08 
in changed the total time less than 4 sec. As already mentioned, 
relatively small changes in descent Mach number produce large changes in 
total time. For descent Mach number changes of ^0.02, time changes greater 
than 12 sec were calculated, and the variation is nearly linear over this 
range. 

4.1.2 Span of Control . - In the previous example, one combination of 
M^/CASjj was determined which would allow the aircraft to fly from the entry 
fix to the metering fix for the specified end conditions. Actually, there 
are an infinite number of combinations of M^j/CAS^j which will accomplish 
this objective. This can be observed in Fig. 4-6 which shows the span of 
control for the B-737. In this figure, the speeds and altitudes given 
earlier were used for the entry and metering fixes and the total distance 
AStot fixed at 76.05 n.mi. The two-gradient wind profile shown on 
Fig. 4-7 for the Denver airport was used and the ground heading was 
0 = 238°. Note that the wind angle given in Fig. 4-7 is 180° from the 
wind direction, and for this example it gives the aircraft side- and tail- 
wind components. Results are given for Segments 4 and 5 the same as 
described earlier, and with these two segments reversed. The descent Mach, 
M(j, was varied from the minimum to the maximum value, and the descent 
calibrated airspeed, CAS(j, was varied over the allowable range. It was 
found that some- combinations of M^/CAS^ produced cross-over altitudes, 
hj^, below hj^j? or greater than h^,. When this occured, the computer 
algorithm changes CAS^ to make h^ = h^^^ when h^ < h^^^ or make 
hx = he when h/ > he* These changes are shown in Fig. 4-8. 

Fig. 4-6 shows that a greater span of control is available when 
Segment 4 is the constant speed segment and Segment 5 is the deceleration 
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seQiTisnt. Rsfs. 2 and 3 had thssG two SGgmsnts rsvsrssd. ThG tVAD snvGlopGS 
merge when because the deceleration segment has zero length. For a 

specified value of At^Q^, the envelopes give the range of which will 
produce the desired profile descent. The corresponding range of CAS^j is 
shown on Fig. 4-9 for Segment 4 the deceleration segment and Segment 5 at 
constant speed. 

Fig. 4-10 gives the envelopes of control spans for aS^q^ = 50, 75, and 

100 n.mi. with no wind. The envelopes for 75 and 100 n.mi. are similar, 

but displaced. However, the envelope for 50 n.mi. is much smaller because 
when < 0.77 a profile descent is not possible for the prescribed end 
conditions. 

4 . . 

4.1.3 Effects of Mind and Non-Standard Day . - Wind can have a signi- 
ficant effect on the time required to perform a profile descent because it 
changes the ground speed. Consider the wind speed (V^) model given in 
Fig. 4-7 for the Denver airport. Calculations were performed for the B-737 
using this wind speed, but the wind direction was for a head wind, side 
wind, and tail wind on a standard day. The envelopes of spans of control 
are shown on Fig. 4-11 for a distance between the entry fix and metering 

fix of 75 n.mi. The influence of the wind is even greater when this dis- 

tance is increased. Note that with the approximations used for Eq. 22, the 
side wind has no effect on the ground speed; thus the calculated profile 
descent is the same as that with no wind. 

It was found that non-standard atmospheric temperatures have a signi- 
ficant effect on calculated profile descents. For the B-737 aircraft 
approaching the Denver airport, profile descents were calculated for sea- 
level temperatures of 59® F (standard day) and 80° F (hot day). 4gain, the 
distance between the entry and metering fixes was 76 n.mi., a required time 
of 702 sec and a descent Mach number of Mj = 0.67 were specified for each 
case. In order to satisfy the end conditions, the CAS^ had to be 290 kts 
for the standard day and 267 kts for the hot day. The results are given 
in Table 4-6. 
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Table 4-6. 

Effect of hot day on 

profile descent. 



STANDARD DAY 

HOT 

DAY 

SEGMENT 

NUMBER 

At^ (sec) 

AS. (n.mi.) 

At^ (sec) 

aS^ (n.mi 

1 

39.7 

3.70 

13.2 

1.18 

2 

89.7 

9.28 

278.0 

27.38 

3 

271.0 

28.03 

140.7 

14.79 

4 

54.8 

5.91 

55.0 

6.20 

5 

244.9 

29.10 

216.9 

26.49 

TOTAL 

700.1 

76.02 

703.8 

76.04 


Mote the differences in the individual segments. It is also of interest 
to note that the fuel consumed for the standard-day profile descent was 
499 lb, whereas it v/as 467 lb for the hot day. These results show that the 
non-standard atmospheric effects are important and should be included in 
the calculations. 

4.1.4 Optimal MM/CASf< Combination for Minimum Fuel . - For fixed 
end conditions, there are normally an infinite number of M^/CAS^j combina- 
tions which could be used to fly the profile descent between the entry and 
metering fixes. Ref. 7 discusses optimal values of CAS^ for minimum fuel 
consumed when the entire descent is performed at constant calibrated air- 
speed. An investigation was made here to determine the optimal combination 
of Mj/CASjj which requires minimum fuel. The computer algorithm calcu- 
lates the fuel consumed for each profile descent possible for increments in 

M. of 0.01 for M . < M , < M . For each value of there is only one 

d mm — d — max 

value of CAS(j that will satisfy the required end conditions. Then the 
computer program determines the M^/CAS^j combination which requires minimum 
fuel consumption. Results for the B-737 flying the profile descent des- 
cribed above with = 75 n.mi. and no wind are shown on Fig. 4-12. 

This figure shows that the differences in total fuel consumption are small. 
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Fuel Weight, 



Figure 4-12. Effect of on fuel consumed. 
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and hence some other criteria would probably be more important in the de- 
termination of the M(j/CAS(j combination. An alternative criterion used in 
the computer program is the determination of the combination of M^/CASd 
for which the difference between the desired metering-fix arrival time and 
minimum arrival time is 30% of the span of time control at that value of 
Mj. The 30% proportionality constant was chosen, based on the assumptions 
that delay, rather than advancement, is prevalent. This scheme allows some 
flexibility to compensate for deviations from the predicted profile and/or 
subsequent changes in the desired metering-fix arrival time generated by 
terminal dynamics. Obviously, this flexibility decreases as the aircraft 
approaches the metering fix. 

4.1.5 Delay Capabilities . - Due to constantly changing conditions 
near an airport, it is necessary to know how much time a flight might be 
delayed at some point in the profile descent and still fly a clean configu- 
ration/idle thrust trajectory to the metering fix. Calculations were per- 
formed for the B-737 approaching the Denver airport to determine the maxi- 
mum delay capability at every position in the profile descent. At any 
position before or in the profile descent this requires changes in the re- 
maining segments that will consume the maximum time to arrive at the mete- 
ring fix at the prescribed altitude and calibrated airspeed there. The 
difference between the maximum time and the time required for the original 
profile descent is the maximum delay capability at that position. 

Fig. 4-13 shows the maximum delay capability as a function of time from the 
metering fix. Prior to arriving at the entry fix, a delay could be 
achieved by initiating the profile descent before the original entry fix. 
For the profile descent planned originally in Fig. 4-13, a maximum delay of 
133 sec is possible when 800 sec from the metering fix. While in Segment 
5, a maximum delay of 103 sec is available and this value does not change 
until Segment 4 begins because the deceleration could be continued beyond 
the planned end-of Segment 5. In Segment 4, the constant speed segment, 
the maximum delay capability decreases rapidly to virtually none at the end 
of Segment 4. 
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It was found that the maximum delay capability available after the 
descent was begun (after Segment 4) was very small. This fact, coupled 
with the additional work load on the pilot to make changes during the 
descent, led to the decision to restrict profile changes to those segments 
before the descent segments. 

4.1.6 Comparison v/ith Flight Simulator . - A profile descent was flown 
on the NASA/LRC B-737 TCV simulator and the results were compared with the 
computational results from the computer algorithm described herein. In 
order to compare with the simulator time and distance for each segment, the 
profile descent was calculated from the conditions listed below: 

(1) Cruise conditions of W = 85,000 lb, h = 34,000 ft, M = 0.78 

(2) Descend with Mj/CASd =0.78/280 kts to h = 10,000 ft 

(3) Decelerate from CAS = 280 kts to 250 kts at 10,000 ft 

(4) Descend to 6,000 ft at 250 kts 

(5) Decelerate from 250 kts to 210 kts at 6,000 ft 

The idle thrust used in the simulator was different from that shown in 
Fig. 4-2. Therefore, the idle thrust used in the calculations was changed 
to duplicate that in the simulator. The results are shown in the table 
below. 

Table 4-7. Comparison of calculated and flight simulator profile descent. 


SEGMENT 

At^ . (sec) 

Ate, (sec) 


ASca, (ft) 

Descent at = 0.78 

25 

24.6 

18,760 

18,814 

Descent at CAS . = 280 kts 
d 

587 

598.6 

379,368 

387,145 

Deceleration at 10,000 ft 

27 

26.3 

13,772 

13,535 

Descent at CAS f 250 kts 

136 

148.8 

64,646 

70,466 

Deceleration at 6,000 ft 

51 

39.9 

21,273 

16,857 

TOTAL 

826 

838.2 

497,819 

506,817 
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The calculated total time and distance was Within 2% of the simulator re- 
sults. An investigation of each segment shows that the calculated results 
for the descent at M^j = 0.78 and deceleration at 10,000 ft segments compare 
well with the corresponding results from the simulator. However, the 
descent segments at constant CAS and the deceleration segment at 6,000 ft 
did not compare nearly so well with the simulator results. Part of the 
differences noted for the two descent segments could be attributed to the 
transition at the beginning and end of these segments. The computer algo- 
rithm assumes instant transition. The larger part of the differences 
between calculated results and simulator results can be traced to pilot 
technique which produces speeds and flight path angles different from the 
prescribed values. 

Of particular note is the effect of flight path angle on the decelera- 
tion segment at 6,000 ft. The simulator results show y = 0.44° at the 
beginning and y = -1.57° at the end of this segment, whereas the desired 
path is Y = 0. Eq. 1 gives the deceleration as 

^ = W -9 ^ • 


Normally the last term on the right would be small compared to the first 
term for small values of Y. However, the first term is small here because 
of idle thrust and the relatively low speed (250 kts CAS). The table below 
shows the magnitude of these two terms for the beginning and end of this 
segment. 

Table 4-8. Comparison of values for deceleration. 


POSITION 

Y 

q(F-D)/W 

-q sin Y 

dV/dt 

beginning 

_0.44° 

-1.22 kt/sec 

-0.15 kt/sec 

-1.37 kt/sec 

end 

-1.57° 

-1.08 kt/sec 

+0.52 kt/sec 

-0.56 kt/sec 
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Note that the flight path angle at the end of this segment produces an 
acceleration (0.52 kt/sec) which is nearly the same magnitude as the 
deceleration (-0.56 kt/sec). If the results above are corrected for Y = 0, 
the simulator results for time and distance in this last segment are very 
close to the value calculated here. Depending on pilot technique, it is 
possible for the regions of positive Y to compensate for the regions of 
negative Y over a deceleration segment; thus producing the same effect as 
Y = 0 throughout the segment. 

4.2 Profile Descent for Heavy -Cl ass of Aircraft 

Previous results were for the B-737 airplane. In this section a pro- 
file descent is described for an aircraft whose characteristics are similar 
to the B-747 aircraft. The pertinent input data to the computer program 
are given below: 

Table 4-9. Input data for heavy-class aircraft profile descent. 


w 

= 

500,000 lb 

s 

= 

5,500 ft^ 

TO 

= 

519° R 

POl 

= 

2116.2 Ib/ft 

HC 

= 

40,000 ft 

HFMAX 

= 

23,000 ft 

HFMIN 

= 

19,000 ft 

HAP 

= 

9,000 ft 

DAP 

= 

36.7 n.mi. 

DSC 

= 

21.2 n.mi. 

VCAP 

= 

170 kts 

VCSC 

= 

210 kts 

VCMF 

= 

250 kts 

AMC 

= 

0.85 

DSTOT 

= 

76.049 n.mi. 

TREQ 

= 

675 sec 

KSEQ = 
NO WIND 

0 
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For this case, the distance from the entry fix to the metering fix is about 
the same as that for one of the B-737 cases (76.049 n.mi.). The cruise 
conditions at the entry fix are = 0.85 and h^ = 40,000 ft. A profile 
descent was required to fly the aircraft from the entry fix to the metering 
fix, where the altitude window is 19,000 ft to 23,000 ft and the calibrated 
airspeed is 250 kts, and arrive there in 675 sec. The aim point is 36.7 
n.irri. from the metering fix at an altitude of 9,000 ft and the speed there 
is 170 kts. A speed control segment is allowed with a calibrated airspeed 
of 210 kts. 

From the conditions specified at the metering fix and aim point, the 
computer algorithm calculated a metering-fix altitude of 22,753 ft. Then 
M^/CASjj = 0.69/296 kts was determined as one combination v/hich would allow 
the aircraft to arrive at the metering fix at the prescribed time of 675 
sec. The additional time required to fly from the metering fix to the aim 
point was 526 sec. 

Results for this profile descent are shown graphically in Fig. 4-14. 
This figure shows the rate of descent, descent angle, and altitude as a 
function of time measured from the entry fix. The rate of descent varies 
nearly linearly with time (and hence altitude) in those segments where des- 
cent is performed at constant calibrated airspeed. However, the rate of 
descent did not change monotonically with time or altitude for Segment 3 
(constant Mach descent). These characteristics for the rate of descent 
were found to hold for other aircraft as well. Therefore, the approximate 
relations given earlier by Eqs. 43 and 45 are justified. It is also ob- 
served that I Y I < 6° and | dY/dt | is small which verifies the approximations 
used in Eq. 9. Although y and dY/dt are discontinuous at the end points of 
descent segments, it was observed in the comparison with the flight simula- 
tor that the effect of these approximations on the arrival times at the 
metering fix and aim point are small. Fig. 4-14 also indicates that flight 
times for a larger aircraft like the B-747 are comparable to those of the 
B-737 for profi4e descents. 
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4.3 Results Using Approximate Relations 


To assess the accuracy of the approximate relations given in Eqs. 37 
through 46 to calculate At^ and aS^ , profile descents were calculated and 
the results are compared below with those from the computer algorithm. The 
results presented here are for the portion of the flight from the entry fix 
to the metering fix. 


EXAMPLE 1: B-737 - h^. = 35,000 ft, M^, = 0.78, = 19,500 ft, 

CAS^r, = 250 kts, AS-rnr = 76.049 n.mi.. At = 650 sec, 

ap TUI req 

e = 238°, winds in Fig. 4-7, M^j/CASj = 0.76/324 kts, 

Tq = 540° R, Po = 2116.2 Ib/ft^ 


Table 4-10. Comparison of calculated and approximate results for B-737. 



At. 

(sec) 

AS. 

(n.mi .) 

SEGMENT 

COMPUTER 

APPROX. 

COMPUTER 

APPROX, 

5 

8.63 

8.38 

1.04 

1.00 

4 

303.28 

307.51 

36.04 

36.24 

3 

155.51 

156.45 

18.80 

18.80 

2 

112.93 

112.65 

13.36 

13.27 

1 

68.50 

67.41 

6.82 

6.74 

TOTAL 

648.85 

652.40 

76.06 

76.05 
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EXAMPLE 2: B-747 - Same case as listed under Table 4-9. 


Table 4-11. Comparison of calculated and approximate results for heavy- 
class aircraft. 


SEGMENT 

Ati 

COMPUTER 

(sec) 

APPROX. 

AS. 

1 

COMPUTER 

(n.mi .) 
APPROX. 

5 

81.4 

79.4 

10.01 

9.74 

4 

130.9 

134.6 

14.39 

14.80 

3 

383.3 

381.7 

43.02 

42.93 

2 

25.9 

25.9 

2.98 

2.98 

1 

53.3 

52.9 

5.64- 

5.59 

TOTAL 

674.8 

674.5 

76.04 

76.04 


Note that the total distance is the input value in each case because AS 4 
is calculated such that the total distance is always correct. The results 
above show that the approximate equations compare well with the computer 
calculations. However, the rates of descent and decelerations at the end 
points of each segment were calculated from the same equations as those 
used in the computer program. Only the region between the end points was 
approximated. When approximate expressions for the rate of descent, 
similar to those used in Refs. 2 and 3, were used significant errors would 
result in some cases. For the B-737 profile descent described in Fig. 15 
of Ref. 2, h^ = 35,000 ft. M^ = 0.78, = 0.62, CAS^ = 250 kts = CAS^^, 

and h|^^ = 19,500 ft. The rate of descent and time for Segment 3 
(Md=0.62) calculated using the data of Ref. 1 are compared with results 
using the approximate relations given by Ref. 2: 
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CALCULATED FROM 

CALCULATED FROM APPROXIMATE EQUATIONS 

PARAMETER DATA IN REF. 1 IN REF. 2 

rate of descent (avg) -2,185 ft/min -2,638 ft/min 

At 3 257 sec 213 sec 

The time required for Segment 3 is in error by 44 sec out of 257 sec. 

Upon reviewing Figure 4 in Ref. 2, it was observed that M(j = 0.62 is 
near the lower end of the curve where the curve-fit could be in error. 

When the approximate equations of Ref. 2 were applied to the profile des- 
cent described above in Example 1, the results compared favorably with the 
present results. For Segment 3 (Mj = 0.76) in that example, a comparison 
of the two methods is given below: 

CALCULATED FROM 

CALCULATED FROM APPROXIMATE EQUATIONS • 

PARAMETER DATA IN REF. 1 IN REF. 2 

rate of descent at be- -3341 ft/min -3329 ft/min 

ginning of Segment 3 

rate of descent at -4320 ft/min -4523 ft/min 

end of Segment 3 

At 3 156 sec 150 sec 

This favorable comparison is obtained because the curve-fit parameter in 
Figure 4 of Ref. 2 is close to calculated values for Mj = 0.76, whereas 
the inaccurate value was found for M^ = 0.62 in that figure. 
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5.0 CONCLUSIONS 


The following conclusions are drawn from the results obtained in this 
report: 

(1) Uncertainty in descent Mach number had the largest influence on 
the time required to fly a profile descent, whereas weight had an 
insignificant influence. 

(2) A greater span of control is available when the initial decelera- 
tion at cruise altitude occurs before the constant speed segment 
as opposed to reversing these two segments. 

(3) Wind and non-standard atmospheric properties have a large effect 
on the time involved in a profile descent and therefore should be 
included in the calculations. 

(4) The optimal combination of descent Mach and calibrated airspeed 
for minimum fuel consumed produces only small savings in fuel. A 
combination based on flexibility for scheduling seems preferable. 

(5) Significant delay capability may be available before the initial 
descent begins, whereas very little is available while in the 
descent. Considering the extra work load on the flight crew, it 
is recommended that profile change commands be given before the 
descent begins. 

(6) Profile descent parameters for heavier aircraft (like the B-747) 
are similar to those of the B-737. 
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(7) Errors from using the point-mass trajectory equations and 
neglecting the transition at the ends of a segment were found to 
be small. Calculated results compared reasonably well with those 
from a B-737 flight simulator. 

(8) The calculations for time and distance in a segment can be sim- 
plified by using approximate relations for the integrals invol- 
ving the rates of descent and deceleration. However, curve fits 
to the rates of descent and deceleration were found to be inaccu- 
rate in some cases. These rates can be calculated accurately 
from the point-mass equations. 
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APPENDIX A 

Pressure and Temperature Variation with Altitude 

In the stratosphere (below 36,089 ft) the standard temperature is 
taken to be (Ref. 4): 

T = Tq - ah (A-1) 

where is sea-level temperature and a is the temperature lapse rate. 

On a standard day 

Tq = 519°R (288°K) 
and 

a = 3.56616 x 10"^ °R/ft 


The pressure is related to the temperature, and hence altitude, by the 
rel ation 


£_ = 
Po 



^5.256114 



5.256114 


(A-2) 


where Pq is sea-level pressure. 

As noted by Knox and Cannon (Ref. 2), nonstandard temperatures and 
pressure affect several Mach numbers, airspeed, and altitude relations used 
in this report. It is assumed that the temperature lapse rate, a, used in 
Eq. A-1 is not affected by nonstandard days and that only the sea-level 
temperature, Tq, changes. Eq. A-2 still holds except Pq must be replaced 
by Pq, the nonstandard sea-level pressure. 

For altitudes greater than 36,089 ft, the standard temperature remains 
constant at Tj^ (temperature at 36,089 ft). In this region, the standard 
pressure is given by Ref. 4 as 


^ = exp [-g(h-36,089)/RT,] (A-3) 

Pi i 


where pj is the pressure given by Eq. A-2 at 36,089 ft altitude. 
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APPENDIX B 

Mach Number, True, and Calibrated Airspeed Relations 


These relations were obtained from Ref. 4. The Mach number is defined 


by 


M = 1 


(B-1) 


where a is the local speed of sound. 




a = W YRT . 


(B-2) 


Here y is the ratio of specific heats (1.4 for air), and R is the gas con- 
stant for air, 

R = 1716 ft^/sec^ / °R 

At subsonic speeds, the ratio of Pitot pressure (p^) to static pres- 
sure (p) is related to the Mach number by 




(B-3) 


This equation can be rearranged to yield 


M‘ 


2 


where 


Pt - P 


Y-1 


Pt - P \ Y 
-rr— + 1 -1 


(B-4) 


„ +1 was substituted for — because airspeed indicators nor- 

P P 


mally measure (p^ - p) instead of p^ alone. Substitute M from Eq. B-1 
into Eq. B-4 to obtain the true airspeed squared as 


2a^ 

Y-1 


Pt - P IT 
^ +1 -1 


(B-5) 
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o 

Since typical airspeed indicators do not measure a and p individually, 
they are calibrated using sea-level values for these parameters. With 
these replacements, the calibrated airspeed squared is obtained from 
Eq. B-5 as 



(B-6) 


The calibrated airspeed differs from true airspeed everywhere except at sea 
level. 

True airspeed can be calculated from calibrated airspeed by substitu- 
ting for (p^ - p) from Eq. B-6 into Eq. B-5 to get 


2 a'" 

w^ry 


P P 

Divide this equation by a^ to obtain 


1 . 


Y 


-1 




2 

Y-1 


1 . 

^0 


Y 


(B-7) 


(B-8) 


relate Mach number, true, and calibrated airspeed. The 


\ L 

Eqs. B-7 and B-8 
Po 

pressure ratio — is obtained from Eq. A-2 as a function of altitude. Note 
that pQ and ag must be standard sea-level values even though the actual 
sea-level values may be nonstandard. 


Cross-Over Altitude 


The cross-over altitude is the altitude which gives the cross-over 
point from the constant Mach descent in Segment 3 to the constant cali- 
brated airspeed- descent in Segment 2. Substitude from Eq. A-2 into 

Po 

Eq. B-8 and solve for the cross-over altitude as: 
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(B-9) 


For nonstandard days, the term in 
plied by: 


square brackets above must be multi- 



1 

5.256114 


to account for nonstandard sea-level pressure in Eq. A-2. Eq. B-9 is 
restricted to h^^ < 36,089 ft. 

If > 36,089 ft, eq. (A-3) must be used to obtain the cross-over 
altitude as 


h 


X 


36,089 ft 



where 


(B-10) 


P " PO P 


and Eq. B-8 is used to obtain — for prescribed values of M and CAS. 
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APPENDIX C 


Relations for dV/dh in Descent Segments 


The rate of climb is given by Eq. 7 which requires a relation for 
dV/dh. This relation is different for constant M and CAS segments. When 
the calibrated airspeed is maintained constant, differentiate Eq. B-7 with 
respect to altitude, using Eqs. A-1, A-2, and B-2 to obtain: 


dV 

dh 


aV 5.256114 a R ^o 
2T " V p~ 


1 - 




Y-1 


1 

Y 



(C-1) 


When the Mach number is maintained constant; differentiate Eq. B-1, 
using Eqs. B-2 and A-1 to obtain: 


dh " " 2T 


(C-2) 


Eqs. C-1 and C-2 are valid for h < 36,089 ft. 

For h > 36,089 ft, T = Tj^ = constant and Eq. A-3 must be used in 
place of Eq. A-2. For constant calibrated airspeed Eq. C-1 must be 
replaced by 
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^ _ JI3. fo 
dh V p 


!o ^^0 

P P 


1 + 


Y-1 

2 




1 

Y 


When the Mach number is constant Eq. C-2 must be replaced by 


dV 

dh 


0 


(C-3) 


(C-4) 


for h > 36,089 ft. 
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APPENDIX D 

Description of Computer Program 


The input data required to calculate the profile descent are: 


W 

S 

TO 

POl 

HC 

HFMAX 

HFMIN 

HAP 

DAP 

DSC 

VCAP 

VCSC 

VCMF 

AMC 

AMMAX 

AMMIN 

VCMX 

VCMN 

DSTOT 

HDGTRK 

TREQ 

KSEQ 


weight (lb) 

wing reference area (ft^) 
sea-level temperature (°R) 
sea-level pressure (Ib/ft^) 
cruise altitude (ft) 
maximum metering-fix altitude (ft) 
minimum metering-fix altitude (ft) 
aim-point altitude (ft) 

distance from metering fix to aim point (n.mi.) 

distance from metering fix to speed control segment (n.mi.) 

calibrated airspeed at aim point (kts) 

calibrated airspeed for speed control (kts) 

calibrated airspeed at metering fix (kts) 

cruise Mach number 

maximum Mach number for aircraft in descent 

minimum Mach number for aircraft in descent 

maximum calibrated airspeed for aircraft in descent (kts) 

minimum calibrated airspeed for aircraft in descent (kts) 

ground distance from entry fix to metering fix (n.mi.) 

track heading from entry fix to metering fix (deg) 

time required to travel from metering fix to entry fix (sec) 

flag; 0 for profile descent as described in this report, 1 for 
Segments 4 and 5 reversed 
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WOMEGA 

WRHO 

WPSI 

W6AM 

WSPDB 

WDIRB 

GROUND 

WGA2 

WRH02 

WPSI2 

WGAM2 

HWP 


surface wind speed at airport (kts) 
wind speed gradient at airport (kts/ft) 
surface wind direction at airport (deg) 
wind direction gradient at airport (deg/ft) 
wind speed bias 
wind direction bias 

altitude of airport above mean sea level (ft) 

wind speed at cruise altitude (kts) 

wind speed gradient at cruise altitude (kts/ft) 

wind direction at cruise altitude (deg) 

wind direction gradient at cruise altitude (deg/ft) 

altitude separating two-gradient wind profile (ft) 


Drag Data: Drag Coefficient CD as a function of Mach number (Afi) 

and lift coefficient (CL). 

Idle-Thrust Data: Engine Thrust (THRUST) as a function of altitude (H) and 

Mach number (AM). 

Fuel Flow Rate: Flow Rate (Ib/sec) of fuel as a function of altitude 

(H), Mach number (AM) and thrust. This information is 
required only if fuel consumption is desired in the cal- 
culation. 


The flowchart for the calculations is given on Figure D-1, and the 
data printed from the computer program are: 


1. Input data 


2. DAP (MAX) 
DT(MAX) 
DAP (MIN) 
DT(MIN) 


maximum distance from metering fix to aim point (n.mi.) 
maximum time for above (sec) 
minimum distance (n.mi.) 
minimum time (sec) 
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Figure D-1. Flov/chart for profile descents. 
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3. , At-j for i = lA, 2A, 3A, and 4A 

hgj, = altitude for speed control (ft) 


4. 


For . 

mm — a — max 


At, 


At. 


(sec) 

max mm ' ' 

iterations for M^j/CAS^j computation 


CASj (kts) 

(ft) (ft) 


ASj (n.mi.)» Ati (sec), and Aw^ (lbs) for i = 1, 


5 


5. Influence coefficients 
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